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It is shown that the resonance Kondo tunneling through a double quantum dot (DQD) with even 
occupation and singlet ground state may arise at a strong bias, which compensates the energy of 
singlet /triplet excitation. Using the renormalization group technique we derive scaling equations 
and calculate the differential conductance as a function of an auxiliary dc-bias for parallel DQD. 



PACS numbers: 72.10.-d, 72.10.Fk, 72.15. Qm, 05.10.Cc 



Many fascinating collective effects, which exist in 
strongly correlated electron systems (metallic compounds 
containing transition and rare-earth elements) may be 
observed also in artificial nanosize devices (quantum 
wells, quantum dots, etc). Moreover, fabricated nanoob- 
jects provide unique possibility to create such conditions 
for observation of many-particle phenomena, which by no 
means may be reached in "natural" conditions. Kondo 
effect (KE) is one of such phenomena. It was found the- 
oretically [jij and observed experimentally [|| that the 
charge-spin separation in low-energy excitation spectrum 
of quantum dots under strong Coulomb blockade mani- 
fests itself as a resonance Kondo-type tunneling through 
a dot with odd electron occupation Af (one unpaired spin 
1/2). This resonance tunneling through a quantum dot 
connecting two metallic reservoirs (leads) is an analog 
of resonance spin scattering in metals with magnetic im- 
purities. A Kondo-type tunneling may be observed un- 
der conditions which do not exist in conventional metal- 
lic compounds. In particular, Kondo effect survives in 
essentially non-equilibrium state when the strong bias 
eV ^> Tk is applied between the leads [|| (Tk is the 
Kondo temperature which determines the energy scale of 
low-energy spin excitations in a quantum dot). The KE 
may be observed as a dynamical phenomenon in strong 
time dependent electric field Q, it may arise at finite 
frequency under light illumination Q . Even the net zero 
spin of isolated quantum dot (even Af) is not the ob- 
stacle for the resonance Kondo tunneling. In this case 
it may be observed in specific types of double quantum 
dots (DQD) [|| or induced by strong magnetic field 0] 
whereas in conventional metals magnetic field only sup- 
presses the Kondo scattering. The latter effect was also 
observed experimentally |q| . 

As was noticed in MQ, quantum dots with even Af 
possess the dynamical symmetry SO (4) of spin rotator 
in the Kondo tunneling regime, provided the low-energy 
part of its spectrum is formed by a singlet-triplet (ST) 
pair, and all other excitations are separated from the ST 
manifold by a gap noticeably exceeding the tunneling rate 
7. A DQD with even in a side-bound configuration 
where two wells are coupled by the tunneling v and only 
one of them (say, /) is coupled to metallic leads (L,R) 



is a simplest system satisfying this condition |6|. Such 
system was realized experimentally in Ref. [ fl0| . 

In the present paper one more unusual property of 
DQDs with even Af is revealed. It is shown that in the 
case when the ground state is singlet \S) and the ST 
gap S 3> Tk, a Kondo resonance channel arises under a 
strong bias eV comparable with 5. The channel opens at 
I eV — 5 1 < Tk, and the Kondo temperature is determined 
by the non- diagonal component Jst = (T\J\S) of effec- 
tive exchange induced by the electron tunneling through 
DQD (Fig. lb). 
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FIG. 1. (a) Double quantum dot in a side-bound configu- 
ration (b) co-tunneling processes in biased DQD responsible 
for the resonance Kondo tunneling. 

The basic properties of symmetric DQD occupied by 
even number of electrons Af — 2n under strong Coulomb 
blockade in each well are manifested already in the sim- 
plest case n — 1, which is considered below. Such DQD 
is an artificial analog of a hydrogen molecule H 2 . If the 
inter-well Coulomb blockade Q is strong enough, one has 
Af = ni + n r , ni = n r = 1, the lowest states of DQD 
are singlet and triplet and the next levels are separated 
from ST pair by a charge transfer gap ~ Q. We assume 
that both wells are neutral at ni, r — 1. Then the effec- 
tive inter-well exchange responsible for the singlet-triplet 
splitting arises because of tunneling V between two wells, 
J = v 2 IQ = S. It is convenient to write the effective spin 
Hamiltonian of isolated DQD in the form 



H i = E 3 \S){S\ + J2 E T\ T v)(Tv\ 



A=S,Tr; 



AA 



(1) 



where X AA ' = |A)(A' 
operator (see, e,g, |ll 



is a Hubbard configuration change 
), E T — E s + 6, r) = ±, are three 



1 



projections of S = 1 vector. Two other terms completing 
the Anderson Hamiltonian, which describes the system 
shown in Fig. la, are 



H b + H t = 



AA ko 



(2) 



The first term describes metallic electrons in the leads 
and the second one stands for tunneling between the leads 
and the DQD. Here a — L,R marks electrons in the left 
and right lead, respectively, the bias eV is applied to 
the left lead, so that the chemical potentials are [Ifl = 
fJ-FR + eV, W is the tunneling amplitude for the well 
I, | A) are one-electron states of DQD, which arises after 
escape of an electron with spin projection er from DQD 
in a state |A). 

In case of strong Coulomb blockade Q ^S> 7, the non- 
equilibrium repopulation of DQD is a weak effect and one 
may start with a second order perturbation calculation in 
tunneling amplitude known as the Schrieffer- Wolff (SW) 
transformation fll}] , where both leads are considered as 
independent subsystems. As is shown in Refs. ||[| the 
SW transformation being applied to a spin rotator results 
in the following effective spin Hamiltonian 



Hint = Et( J ^ S + &) ■ W) + J%X 

( \ n 
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(3) 



Here s aa >— J2kk' c kaa^ Ck ' a 
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f, 1 are the Pauli matrices and unity matrix respec- 
tively. The effective exchange constants are 



J, 



AA' 



w 
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2 \e Fa -(E s /2 + S) e Fa , - (E 8 /2 + 5) 
Two vectors S and P with spherical components 

S+ = V2 (A 10 + A " 1 ) , S~ = V2 (A 01 + X- 10 ) , 
5 Z = A 11 -A- 1 - 1 , P z = - (X os + X so ) , 
P+ = V2 (A 1S - A 5 " 1 ) ,P~ = V2 (A S1 - A~ 1S ) . (4) 

obey the commutation relations of 04 algebra 

[Sj,Sk] = iejkiSi, [Pj,Pk] = iejkiSi, [Pj,S k ] = ie jk iPi 

(J, k, I are Cartesian coordinates, ejki is a Levi-Chivita 
tensor). These vectors are orthogonal, S • P = 0, and 
the Casimir operator is S 2 + P 2 = 3. Thus, the singlet 
state is involved in spin scattering via the components of 
the vector P. 

We use SU(2)-\ike semi-fermionic representation for S 
operators @,[L3| 



S+ = V2(/t/-i + flfo), S- = V2(f f _J Q + flh), 

Sz = fth-fl 1 f-i, (5) 



where /J. are creation operators for fermions with spin 
"up" and "down" respectively, whereas /o stands for 
spinless fermion |l^,^3| . This representation can be gen- 
eralized for 5*0(4) group by introducing another spinless 
fermion f s to take into consideration the singlet state. 
As a result, the P-operators are given by the following 
equations: 

P+ = V2(flf s - /t/_!), P- = V2(/t/i - flj a ), 



P* = -Wf. + ftfo). 



(6) 



The Casimir operator S 2 +P 2 = 3 transforms to the local 
constraint £a=± A s /a/a = 1. 

The spin Hamiltonian is now given by 
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where S d and P d (d=x,y,z) are 4x4 matrices defined by 
relations (§ - §) and J s = J S5 , J T = J TT and J ST are 
singlet, triplet and singlet-triplet coupling SW constants, 
respectively. 

To develop the perturbative approach for T > Tk 
we introduce the temperature Green's functions (GF) 
for electrons in a dot, Ga(t) = — (T t /a(t)/ a (0)), and 
GF of left (L) and right (R) electrons in the leads 
G L ,R{k,r) = -(T T c LMa (k,T)c i L Ra (k,0)). Performing a 
Fourier transformation in imaginary time for bare GF's, 
we come to following expressions: 

G L( e ") = ( ie n ~ £fe +A*i,fi) _1 , 



g° s (e n ) = {ie n -E s )- x , 



1,0,1 



(8) 



with e„ = 2ttT(?i + 1/2) and u) m = 2i:T{m+ 1/3) 

The first leading and next to leading parquet diagrams 

are shown on Fig. 2. 

In equilibrium state eV = the elastic Kondo tun- 
neling arises only provided Tk 3> S in accordance with 
the theory of two- impurity Kondo effect dJl^l . Now we 
will show that in the opposite limit Tk <C <5 the elastic 
channel emerges at eV ~ 5. 

Corrections to the singlet vertex T(u>, 0; oj, eV) are cal- 
culated using an analytical continuation of GF's to the 
real axis u> and taking into account the shift of the chem- 
ical potential in the left lead. In a weak coupling regime 
T > Tk the leading non-Born contributions to the tunnel 
current are determined by the diagrams of Fig. 2 b-e. 



2 



\ y \ / \ y 





d) 



5) \ 




FIG. 2. Leading (b,d) and next to leading (c,e) parquet 
diagrams determining renormalization of J s (a). Solid lines 
denote electrons in the leads. Dashed lines stand for electrons 
in the dot. 

The effective vertex shown in Fig. 2b is given by the 
following equation 



T ST jTS 
J LL J LR 



E 1 



f(e k L - eV) 



(9) 



Changing the variable ekL for ekL — eV one finds that 



-,(26), 



JllJl^rx (D/maK{u,(eV-5),T}). Here 
D ~ ef is a cutoff energy determining effective band- 
width, ^ is a density of states on a Fermi level and 
/(e) is the Fermi function. Therefore, under condition 
\eV — S\ <C maxfeV, 5] this correction does not depend on 
eV and becomes quasielastic. 

Unlike the diagram Fig. 2b, its "parquet counterpart" 
term Fig. 2c contains eV + S in the argument of the 
Kondo logarithm: 



-(2c)/ A _ T. 

lr yu) - J, 



ST T TS 
LL J LR ' 



fjekL ~ eV) 

- tkL + HL + 



(10) 



At eV ~ S ^S> T, uj this contribution is estimated as 



LR. 



JE&te (D/(eV + 6)) « r^( W ). 



Similar estimates for diagrams Fig. 2d and 2e give 



,(2d) 



(W) ~ Jj 



ST T T 



LL u LL 



Jt r Jllv 2 In (fl/max{ W , (eV - 5), T}) 



( W ) ~ J s L T L J T LL Jr R v 2 In (D/max^, (eV - T}) x 
x In (L>/max{w, e y,T}. (11) 

Then r^(w) < r^(w) at eV 5. 

Thus, the Kondo singularity is restored in strongly 
non-equilibrium conditions when the energy loss 5 in a 
singlet-triplet excitation is compensated by the external 
voltage applied to the lead, but the leading sequence of 
most divergent diagrams degenerates in this case from a 
parquet to a ladder series. 

Following the poor man's scaling approach, we derive 
the system of coupled renormalization group (RG) equa- 
tions for (0). The equations for LL co-tunneling are: 



FIG. 3. Irreducible diagrams contributing to RG equations. 
Hatched boxes and circles stand for triplet-triplet and sin- 
glet-triplet vertices respectively. Notations for lines are the 
same as in Fig. 2 
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The scaling equations for j£ R are as follows: 



(12) 
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(13) 



One-loop diagrams corresponding to the poor man's scal- 
ing procedure are shown in Fig. 3. To derive these equa- 
tions we collected only terms ~ ( J T )" ln n+1 (D/T) ne- 
glecting contributions containing ln[D/(eV)]. The anal- 
ysis of RG equations beyond the one loop approximation 
will be published elsewhere. 



The solution of the system (13) reads as follows: 
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(14) 



Here a — L, a' = L, R. The Kondo temperature is 
determined by triplet-triplet processes only. It is given 
by T K = Dexp[-l/{vj£)]. 

The differential conductance G{eV,T)/G^ is the uni- 
versal function of two parameters T/Tx and V/Tx (see 
Fig. 4), Go = e 2 /irh: 



G/G ~ In" 2 (max[( e y - S),T]/T K ) 



(15) 



Finite decoherence rate T1/T4 effects discussed in ]16j in 
a context of strongly-nonequilibrium transport through 
QD with S = 1/2 do not arise in our case. According 
to the Non-Crossing Approximation (NCA) description, 
the origin of Td is inelastic spin relaxation of Kondo state. 
In the problem considered the ground state is S = sin- 
glet and the spin-relaxation is absent. A Kondo-channel 



3 



arises only in virtual states of L-R co-tunneling. Repop- 
ulation effects at eV > 5 should result in asymmetry 
of Kondo-peak [[lTJ, but this effect is beyond our quasi 
equilibrium approach. 




FIG. 4. The Kondo conductance as a function of dc-bias 
eV/T K and T/T K . The singlet-triplet splitting 5/T K = 10. 

Thus, we have shown that the tunneling through sin- 
glet DQDs with 5 3> Tk exhibits a peak in differential 
conductance at eV ~ S instead of the usual zero bias 
Kondo anomaly (see Fig. 4) which arises in the opposite 
limit, 5 < Tk. Therefore, in this case the Kondo effect in 
DQD is induced by a strong external bias. The scaling 
equations (|l3|), ( |l4| ) can also be derived in Schwinger- 
Keldysh formalism (sec |L3| and also jl6|]) by applying 
the "poor man's scaling" approach directly to the dot 
conductance The detailed analysis of the model (Q) 
in a real-time formalism will be presented elsewhere. 

We discuss yet another possible experimental real- 
ization of resonance Kondo tunneling driven by exter- 
nal electric field. Applying the alternate field V = 
V ac cos(wt) to the parallel DQD, one takes into consider- 
ation two effects, namely (i) enhancement of Kondo con- 
ductance by tuning the amplitude of ac- voltage to satisfy 
the condition \eV ac — 8\ -C Tk and (ii) spin decoherence 
effects due to finite decoherence rate Q . One can expect 
that if the decoherence rate h/r ^> Tk, 



G peak /G ~ In" 2 (Tl/tTk) 
whereas in the opposite limit Ti/t <C Tk, 



Gpeak = G{V ac COS [uXt]) 



(17) 



is averaged over a period of variation of ac bias. In this 
case the estimate ( |l5| ) is also valid. 

In conclusion, we have provided the first example of 
Kondo effect, which exists only in non-equilibrium con- 
ditions. It is driven by external electric field in tunneling 
through a quantum dot with even number of electrons, 
when the low-lying states are those of spin rotator. This 
is not too exotic situation because as a rule, a singlet 
ground state implies a triplet excitation. If the ST pair 
is separated by a gap from other excitons, then tuning the 



dc-bias in such a way that applied voltage compensates 
the energy of triplet excitation, one reaches the regime 
of Kondo peak in conductance. This theoretically pre- 
dicted effect can be observed in dc- and ac-biased double 
quantum dots in parallel geometry. 

This work is partially supported (MK) by the Eu- 
ropean Commission under LF project: Access to the 
Weizmann Institute Submicron Center (contract num- 
ber: HPRI-CT-1999-00069). The authors are grateful to 
Y. Avishai, A. Finkel'stein, Y. Gefen, H. Kroha, A. Rosch 
and M.Heiblum for numerous useful discussions. The fi- 
nancial support of the Deutsche Forschungsgemeinschaft 
(SFB 410) is gratefully acknowledged. The work of KK 
is supported by ISF grant. 



[1] 
[2] 
[3] 
[4] 

[5] 



[6] 
[7] 



(16) [8] 



[9] 
[10] 

[11] 

[12] 
[13] 

[14] 
[15] 



L.I.Glazman and M.E.Raikh, JETP Lett. 47, 452, 
(1988); T.K.Ng and P.A.Lee, Phys.Rev.Lett., 61, 1768 
(1988). 

D. Goldhaber-Gordon et.aL, Phys.Rev.Lett. 81, 5225 
(1998); S.M. Cronenwett et al, Science, 281, 540 (1998); 
F. Simmel et al, Phys. Rev. Lett. 83, 804 (1999). 
Y. Meir et al, Phys.Rev.Lett. 68, 2512 (1992), Phys.Rev. 
B49, 11040 (1994); T.K. Ng, Phys. Rev. Lett. 76, 487 
(1996); M.H. Hettler et al, Phys. Rev. B58, 5649 (1998); 
Y. Goldin and Y. Avishai, Phys. Rev. Lett. 81, 5394 
(1998), Phys. Rev. B61, 16750 (2000); A. Kaminski et 
al, Phys. Rev. B62, 8154 (2000); R. Lopez et al, Phys. 
Rev. B 64, 075319 (2001). 

K. Kikoin and Y. Avishai, Phys. Rev. B62, 4647 (2000); 
T.V. Shahbazyan et al, Phys. Rev. Lett., 84, 5896 (2000); 
T. Fujii and N. Kawakami, Phys. Rev. B63, 054414 
K.Kikoin and Y. Avishai, Phys.Rev.Lett. 86, 2090 (2001). 
M. Pustilnik et al, Phys. Rev. Lett. 84, 1756 (2000); M. 
Eto and Yu. Nazarov, Phys. Rev. Lett. 85, 1306 (2000), 
Phys. Rev. B64, 085322 (2001); D. Giuliano et al, Phys. 
Rev. B63, 125318 (2001). M. Pustilnik and L.I. Glaz- 
man, Phys. Rev. Lett. 85, 2993 (2000), Phys. Rev. B64, 
045328 (2001). 

S. Sasaki et al, Nature 405, 764 (2000); J. Nygard et al, 
Nature 408, 342 (2000). 

K.Kikoin and Y. Avishai, Phys.Rev. B65, 115329 (2002). 
L.W.Molenkamp et al, Phys. Rev. Lett. 75, 4282 (1995) 

A. C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, 1993). 
V.Popov and S.Fedotov, Sov.Phys.JETP 67, 535 (1988). 
M.N.Kiselev and R.Oppermann, Phys. Rev. Lett 85 
(2000), M.N.Kiselev et al, Eur. Phys. J B22, 53 (2001). 

B. A.Jones and C.M.Varma, Phys.Rev. B40, 324 (1989). 
Possibility of additional Kondo peaks at eV 3> Tk in 
strong magnetic field was noticed also in the first and 
last papers of Ref . [7] . 



[16] A.Rosch et al, bond- mat /0202404 O.Parcollet and 



C.Hoolcy, |cond-mat/0202425 



cond-mat/0203001, 3ond-mat/0205004 



P.Coleman and W.Mao, 



4 



